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The total B — > X s y decay rate and the CKM-matrix element \V„b\ play an important role in finding 
indirect evidence for new physics affecting the flavor sector of the Standard Model, complemen- 
tary to direct searches at the LHC and Tevatron. Their determination from inclusive B-meson 
decays requires the precise knowledge of the parton distribution function of the b quark in the 
B meson, called the shape function. We implement a new model-independent framework for the 
shape function with reliable uncertainties based on an expansion in a suitable set of basis func- 
tions. We present the current status of a global fit to BaBar and Belle data to extract the shape 
function and the B — > X s y decay rate. 
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Figure 1: On the left-hand side a SM graph for the lowest order B — > X s y decay is shown. In new-physics 
models with two Higgs doublets such as the minimal supersymmetric Standard Model the charged W -boson 
can be replaced by a charged Higgs, as shown on the right-hand side. 



1. Searching for new physics in the flavor sector 

If physics beyond the Standard Model (SM) is present in the flavor sector, it can modify flavor 
changing neutral current processes through the exchange of unknown virtual particles. The inclu- 
sive B — ^ X s y decay via a radiative penguin diagram proves very sensitive to such corrections in 
many new-physics models, which could manifest itself in a modification of the total decay rate. 
The lowest order B — > X s / decay and a new-physics contribution are shown in Fig. 1. To constrain 
such contributions, a good knowledge of the SM prediction and precise measurements are needed. 
Hence, a model-independent framework, as proposed in Ref. [1], that combines all available in- 
formation into a global fit is desirable to accurately determine the B —>X s y decay rate. Here we 
present the current status and preliminary results of such a global fit to the available B — > X s y data. 



2. The B ->X S y decay 



The current analysis of the experimentally measured B — > X s y branching fraction is per- 
formed by extrapolating the measurements to the partial branching fraction for E y > 1.6 GeV, 
which yields [2] &(E y > 1.6) = (3.55 ± 0.24 ± 0.09) x 10~ 4 . This value is then compared to 
the fixed next-to-next-to-leading order (NNLO) SM prediction from Refs. [3, 4], {%(E y > 1.6) = 
(3. 15 ±0.23) x 10~ 4 . The extrapolation still requires a theoretical calculation of the decay rate 
with the measured Ey cut. It assumes a model for the &-quark distribution function, the shape func- 
tion, which introduces an unknown systematic uncertainty. In addition, the extrapolation only uses 
a single branching fraction from each experimental analysis, which is usually the one that has the 
smallest cut on Ey and correspondingly the largest experimental uncertainty. This means that only 
a small subset of all experimental information is used, and more precise measurements at higher 
values of Ey are not utilized. 

A method that avoids these drawbacks was proposed in Ref. [1]. Starting from the effective 
electroweak Hamiltonian, the B — > X s y differential decay rate can be written as 

drv^ 1 = ^ £ 3 (m ,s )2|vrf 

+ YW m F m (m B -2Ey) +£^(cf lcl C i )+ £ 0(QCj). (2.1) 

m J #7 i,#7 

The complete expressions entering Eq. (2.1) will be given in Ref. [5]. Here, F(k) is the nonper- 
turbative contribution to the Z?-quark distribution function in the B meson. Its precise definition is 
given in Ref. [1]. It is F(k) that determines the shape of the B — > X s y spectrum at large Ey and 
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Figure 2: Left: The first five basis functions f„(x). Right: Illustration of the convergence of the expansion. 
The black line is a Gaussian model function F(k) and (k) are its expansions including up to N + 1 terms. 

which we will refer to as the shape function. The function W77 (k) contains the perturbative correc- 
tions to the spectrum resummed to next-to-next-to-leading-logarithmic order, and including the full 
NNLO corrections. At lowest order, Wjj(k) = 8(k). The F m (k) are 1/m^ suppressed subleading 
shape functions. In a fit to B — ^ X s 7 data only, they can be absorbed into F{k) at lowest order in 
a s . The terms proportional to are included at next-to-leading order for i = 1,2,8 using the 
SM values for Ci 2,8- They have almost no effect on the fit, because they are very small in the 
experimentally accessible region of the photon energy spectrum. 

The coefficient C 7 ncl = Cf(no)m b (no) /m\ s + ■■■, where Cf (jUn) is the standard effective Wil- 
son coefficient multiplying O7 = e/(l67i 2 )m b sa^ lv F llv PRb in the electroweak Hamiltonian. The 
ellipses denote all virtual contributions from other operators that generate the same effective b^-sy 
vertex. This includes the terms that cancel the /^j dependence of Cy ff (/Xo)mi,( i Uo), such that C 7 ncl 
is /Xo-independent to the order one is working at. In our analysis we consider |C 7 ncl V^V f *| as the 
unknown parameter that parametrizes the total B — > X s 7 rate. It is extracted simultaneously with 
F(k) from a fit to the measured photon energy spectra. This fitted value can then be compared to 
its SM prediction. In this approach all measurements contribute to constrain the B — > X„ 7 rate. 

3. Extraction of the shape function F(k) from data 

The shape function can be expanded in a complete, orthogonal basis f n {x), constructed in 
Ref. [1] and illustrated in Fig. 2 (where A ~ Aqcd is a parameter of the basis): 
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F(k) 



/ k\ f 
£ Cnfn (y) with / AkF(k) = £ C \ = 1 . 



(3.1) 



The shape of F(k) is then parametrized by the basis coefficients c n , which can be fitted from data. 
The experimental uncertainties and correlations in the measured spectra are captured in the uncer- 
tainties and correlations of the fitted coefficients c„. In practice, the expansion must be truncated 
after ,/V terms, which introduces a residual model dependence from the chosen functional basis. The 
overall size of this truncation uncertainty scales as 1 — Y*n=o c n an ^ tne truncation order ,/V should 
be chosen such that this uncertainty is small in comparison to the experimental uncertainties of the 
fitted coefficients. In other words, we allow the available data to determine the precision to which 
the functional form of the shape function is known, by including as many basis coefficients in the 
fit as are required by the precision of the data. Hence, this approach allows for an experimental 
determination of the shape function with reliable and completely data-driven uncertainties. 
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Figure 3: The Be/fe and BaBar spectra from Refs. [6, 7, 8]. The histograms show the result of the fit. 
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Figure 4: Left: The extracted F(k) (with absorbed 1/m^ corrections) for two (co,i), three (co,i,2), four 
(^0,1,2,3), and five (00,1.2,3,4) coefficients and basis parameter X = 0.5 GeV. The colored envelopes are 
given by the uncertainties and correlations of the extracted coefficients c„. Right: The extracted values 
of |C™ cl VffcVf*! and m l b s with their respective Aj 2 = 1 contours. The grey band shows the SM value for 

\c? cl v tb v*\. 



4. Fit results 

Using the above framework, we perform a global x 2 fit for F(k) and |C" lcl V,£V,*| to all available 
B — > X S J photon energy spectra: The recent Belle measurement in Ref. [6], and the two BaBar 
measurements in Refs. [7, 8]. We have extensively tested our fitting procedure using pseudo- 
experiments. So far, the fit includes experimental uncertainties only. A preliminary study indicates 
that the theoretical uncertainties in the fit results are about the same size as the experimental ones. 

The fit result using four basis coefficients and X = 0.5 GeV as the basis parameter is shown 
in Fig. 3. The fit has a # 2 /dof = 27.67/38 and describes the measured spectra very well. The 
extracted shape function from fits with two to five coefficients are shown on the left in Fig. 4. On 
the right, we show the corresponding results for V^V^I and ni l b s , where the latter is computed 
from the moments of the extracted shape function. The fitted value for |Cy ncl V^V r *| agrees within 
one standard deviation with its SM value, using |V (6 V,* I = 40.68+^ x 1(T 3 and the next-to-leading 
order SM prediction for Ci, ncl = 0.354+q qJ2- For a more stringent comparison an evaluation of Ci, ncl 
in the SM at NNLO using Refs. [3, 4] would be very useful. 
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The results in Fig. 4 show a convergent behavior as the number of basis functions is increased. 
The larger fit uncertainties with more coefficients originate from the larger number of degrees of 
freedom in the fit. A reliable value for the experimental uncertainties is given when the central 
values have converged and the last coefficient, here C4, is compatible with zero within its uncer- 
tainties. At this point the truncation uncertainty can be neglected compared to the fit uncertainties 
in the coefficients. This also implies that using a fixed model function and fitting one or two model 
parameters will in general underestimate the true uncertainties of the shape function. We have also 
checked that using different basis parameters A = 0.4 GeV and 0.6 GeV yields consistent results. 

5. Conclusions 

We presented preliminary results from the first global fit for the nonperturbative shape function 
and the B -^X s y decay rate, parametrized by \Of c{ V t b V t *\, within a model-independent framework. 
The extracted value of \C l f c{ V t b V t *\ agrees with the SM prediction. From our fit we also determine 
m\ s . This constitutes the first step towards a global fit that combines all available data on both 
B X s y and B -> XJvi to determine \Cf cX V tb V*\ and \V ub \ in a model-independent way. 
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